%% Creat  mateix features for relieff
clc
clear all
cd ('E:')
load ('Deltapre_TDCSinput')
load ('Thetapre_TDCSinput')
load ('Deltapre_ECCinput')
load ('Deltapre_Bothinput')
load ('Deltapre_Noneinput')

load ('Thetapre_TDCSinput')
load ('Thetapre_TDCSinput')
load ('Thetapre_ECCinput')
load ('Thetapre_Bothinput')
load ('Thetapre_Noneinput')

load ('Alphapre_TDCSinput')
load ('Alphapre_TDCSinput')
load ('Alphapre_ECCinput')
load ('Alphapre_Bothinput')
load ('Alphapre_Noneinput')

load ('Betapre_TDCSinput')
load ('Betapre_TDCSinput')
load ('Betapre_ECCinput')
load ('Betapre_Bothinput')
load ('Betapre_Noneinput')

load ('Gammapre_TDCSinput')
load ('Gammapre_TDCSinput')
load ('Gammapre_ECCinput')
load ('Gammapre_Bothinput')
load ('Gammapre_Noneinput')
%% Delta
% 
relieffinput=vertcat(Deltapre_TDCSinput,Deltapre_ECCinput,Deltapre_Bothinput,Deltapre_Noneinput);
y=[1;1;1;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;4;4;4];%interventions (1=Before tdcs , 2= before ECC, 3=before Both, 4= before None)
k=10;
[idx,weights] = relieff(relieffinput,y,k)
save idx_Delta idx

%% Theta
relieffinput=vertcat(Thetapre_TDCSinput,Thetapre_ECCinput,Thetapre_Bothinput,Thetapre_Noneinput);
y=[1;1;1;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;4;4;4];%interventions (1=Before tdcs , 2= before ECC, 3=before Both, 4= before None)
k=10;
[idx,weights] = relieff(relieffinput,y,k)
save idx_Theta idx
%% Alpha
relieffinput=vertcat(Alphapre_TDCSinput,Alphapre_ECCinput,Alphapre_Bothinput,Alphapre_Noneinput);
y=[1;1;1;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;4;4;4];%interventions (1=Before tdcs , 2= before ECC, 3=before Both, 4= before None)
k=10;
[idx,weights] = relieff(relieffinput,y,k)
save idx_Alpha idx
%% Beta
relieffinput=vertcat(Betapre_TDCSinput,Betapre_ECCinput,Betapre_Bothinput,Betapre_Noneinput);
y=[1;1;1;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;4;4;4];%interventions (1=Before tdcs , 2= before ECC, 3=before Both, 4= before None)
k=10;
[idx,weights] = relieff(relieffinput,y,k)
save idx_Beta idx
%% Gamma
relieffinput=vertcat(Gammapre_TDCSinput,Gammapre_ECCinput,Gammapre_Bothinput,Gammapre_Noneinput);
y=[1;1;1;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;3;3;3;3;3;3;3;3;3;3;3;3;4;4;4;4;4;4;4;4;4;4;4;4];%interventions (1=Before tdcs , 2= before ECC, 3=before Both, 4= before None)
k=10;
[idx,weights] = relieff(relieffinput,y,k)
save idx_Gamma idx